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Margaret  M.  Murnane1  and  Charles  G.  Durfee14 


We  generate  angularly  isolated  beams  of  circularly  polarized  extreme  ultraviolet  light  through  the  first  implementation  of 
non-collinear  high  harmonic  generation  with  circularly  polarized  driving  lasers.  This  non-collinear  technique  offers 
numerous  advantages  over  previous  methods,  including  the  generation  of  higher  photon  energies,  the  separation  of  the 
harmonics  from  the  pump  beam,  the  production  of  both  left  and  right  circularly  polarized  harmonics  at  the  same 
wavelength  and  the  capability  of  separating  the  harmonics  without  using  a  spectrometer.  To  confirm  the  circular 
polarization  of  the  beams  and  to  demonstrate  the  practicality  of  this  new  light  source,  we  measure  the  magnetic  circular 
dichroism  of  a  20  nm  iron  film.  Furthermore,  we  explain  the  mechanisms  of  non-collinear  high  harmonic  generation  using 
analytical  descriptions  in  both  the  photon  and  wave  models.  Advanced  numerical  simulations  indicate  that  this  non-collinear 
mixing  enables  the  generation  of  isolated  attosecond  pulses  with  circular  polarization. 


The  high  harmonic  generation  (HHG)  process  can  be  used  to 
produce  attosecond  bursts  of  extreme  ultraviolet  (EUV)  and 
soft  X-ray  light  using  tabletop -scale  lasers1-3,  enabling 
numerous  applications,  including  nanoscale  imaging4-8,  studies 
of  nanoscale  energy  transport9  and  element- specific  probes  of 
charge  and  spin  dynamics  on  the  femtosecond  timescale10-14. 
HHG  produces  EUV  light  through  an  ionization-and-recollision 
process15,  and  because  recollision  does  not  occur  when  an  atom 
is  irradiated  with  circularly  polarized  light,  for  many  years  it 
was  believed  to  be  impossible  to  generate  circularly  polarized 
EUV  light  by  the  HHG  process.  This  precluded  many  powerful 
techniques  (such  as  EUV  and  soft  X-ray  magnetic  circular  dichro¬ 
ism  (MCD)  and  photoelectron  circular  dichroism)  from  being 
implemented  on  the  tabletop.  Consequently,  there  have  been 
numerous  experimental16-18  and  theoretical19,20  efforts  to  directly 
produce  circularly  polarized  harmonics  using  a  single  driving 
laser  beam,  but  such  efforts  have  yielded  only  elliptically  polarized 
high  harmonics. 

In  a  recent  breakthrough,  Fleischer  et  al.21  experimentally 
demonstrated  that  HHG  driven  by  two  collinear  beams  with 
counter-rotating  circular  polarization  can  directly  generate  EUV 
light  with  circular  polarization,  as  was  predicted  theoretically22,23. 
Moreover,  phase  matching  of  collinear  circularly  polarized  HHG 
is  remarkably  robust,  providing  a  bright  source  of  EUV  light  and 
thereby  enabling  tabletop  EUV  MCD  measurements24.  However, 
collinear  circularly  polarized  HHG  requires  driving  lasers  at  differ¬ 
ent  wavelengths21,24,25,  which  increases  the  experimental  complexity. 
Furthermore,  collinear  HHG  produces  a  pulse  train  of  three  linearly 
polarized  attosecond  bursts  per  laser  cycle,  where  the  polarization 
rotates  120°  between  each  burst22  (and  sums  to  produce  circular 
polarization  at  the  sample).  Consequently,  if  very  short- duration 
driving  pulses  were  used  to  produce  an  isolated  attosecond  pulse, 
it  would  have  linear  (or  in  some  cases  elliptical26)  polarization. 


In  this  work  we  demonstrate  a  new  approach:  non-collinear  cir¬ 
cularly  polarized  HHG  (NCP-HHG),  in  which  circularly  polarized 
high  harmonics  are  generated  using  crossed  beams  of  counter-rotat¬ 
ing  circularly  polarized  driving  lasers  (Fig.  la).  NCP-HHG  offers 
several  key  benefits  not  found  in  other  HHG  methods 
(Supplementary  Fig.  1).  First,  circularly  polarized  harmonics  can 
be  generated  using  a  single-colour  driving  laser,  maximizing  the 
photon  energies  obtained.  Second,  angularly  isolated  circularly 
polarized  beams  of  both  left  and  right  helicity  can  be  generated  at 
the  same  wavelength,  enabling  precision  differential  measurements 
of  circular  dichroism  at  a  single  wavelength  in  a  single  exposure. 
Third,  by  adjusting  the  crossing  angle  of  the  driving  laser  beams, 
the  different  harmonic  orders  can  be  angularly  separated  without 
the  need  for  a  spectrometer  (Fig.  la,  Supplementary  Fig.  2). 
Finally,  the  harmonics  are  naturally  separated  from  the  driving 
laser  beams,  allowing  a  fragile  sample  to  be  placed  directly  into 
the  EUV  beams  without  the  need  for  filters,  thus  significantly 
increasing  the  usable  HHG  flux.  Here,  we  use  the  NCP-HHG 
source  to  perform  an  EUV  MCD  measurement  of  a  20-nm-thick 
iron  sample,  confirming  that  the  EUV  beams  are  circularly  polar¬ 
ized  with  opposite  helicity  and  demonstrating  that  the  brightness 
and  stability  of  the  source  are  sufficient  to  measure  the  ultrafast 
dynamics  of  magnetic  materials. 

We  explain  the  mechanisms  of  non-collinear  HHG  using  both 
intuitive  physical  models  as  well  as  advanced  numerical 
calculations.  In  the  photon  picture  (Fig.  lb),  we  show  that  the 
NCP-HHG  method  is  governed  by  simple  selection  rules,  allowing 
the  emission  angles  of  EUV  emission  to  be  calculated  analytically. 
Using  a  simple  wave-mixing  model  (Fig.  lc),  we  find  that  the 
circularly  polarized  high  harmonics  are  generated  by  driving 
fields  that  are,  in  fact,  locally  linear.  The  linearity  of  the  field 
experienced  by  each  atom  enables  the  optimal  cutoff  photon 
energies  of  linearly  polarized  HHG  to  be  realized  for  circularly 
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Photon  model:  addition  of  photon  wavevectors 


Figure  1  |  Non-collinear  circularly  polarized  high  harmonic  generation  (NCP-HHG).  a,  In  the  NCP-HHG  experiment,  two  counter-rotating  circularly  polarized 
femtosecond  laser  pulses  are  focused  into  a  noble  gas  (Xe,  Ar  or  Ne)  to  produce  both  left  and  right  circularly  polarized  beams  of  EUV  light  via  HHG.  The 
detector  shows  experimental  data  for  NCP-HHG  driven  by  267  nm  lasers,  where  EUV  beams  at  23  eV  (inner)  and  14  eV  (outer)  emerge  at  different  angles,  well 
separated  from  the  driving  laser  beams,  b,  In  the  photon  model,  the  output  direction  of  each  harmonic  follows  the  simple  vector  addition  of  the  wavevectors  for 
each  absorbed  photon,  c,  In  the  wave  model,  the  two  circularly  polarized  beams  sum  to  yield  an  electric  field  that  exhibits  linear  polarization,  which  rotates 
as  a  function  of  the  transverse  position  across  the  laser  focus,  producing  a  'rotating  polarization  grating'  that  generates  circular  polarization  in  the  far  field. 
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Figure  2  |  NCP-HHG  using  two  800  nm  pulses,  a,  The  EUV  output  of  the  NCP-HHG  process  is  recorded  after  passing  through  a  spectrometer,  which 
disperses  the  harmonics  in  the  vertical  direction.  Each  harmonic  consists  of  two  beams,  one  at  an  angle  left  of  the  centreline  (left  circular  polarization)  and 
one  to  the  right  of  the  centreline  (right  circular  polarization).  b,c,  Numerical  HHG  simulations  driven  by  4.3-cycle  800  nm  laser  pulses  demonstrate  that 
the  left  beamlets  are  left  circularly  polarized  and  the  right  beamlets  are  right  circularly  polarized,  d,  The  divergence  of  each  harmonic  order  q  follows 
conservation  of  momentum,  which  predicts  that  the  angle  of  harmonic  emission,  0q,  is  related  to  the  crossing  angle  of  the  two  laser  beams,  0 1f  by 
sin  6q  =  ±sin  0^/q.  e,  The  temporal  structure  of  the  harmonic  emission  is  a  pulse  train  of  many  circularly  polarized  attosecond  bursts. 
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Figure  3  |  Photon  and  wave  models  of  non-collinear  HHG.  a,  For  the  mixing  of  linearly  polarized  laser  beams,  the  only  constraint  is  that  the  total  number  of 
absorbed  photons  is  odd.  This  allows  many  combinations  of  photons  for  any  given  harmonic  (in  this  case  the  seventh  harmonic)  and  results  in  EUV  beams 
emerging  at  many  different  angles,  b,  In  the  wave  model,  two-beam  mixing  produces  an  amplitude-modulated  grating  at  the  focus,  which  produces  multiple 
peaks  in  the  far  field,  c,  For  counter-rotating  circularly  polarized  beams,  conservation  of  spin  angular  momentum  requires  that  nR  =  nL±  1,  restricting  the 
output  to  only  two  beams  for  each  harmonic,  d,  In  the  wave  picture,  the  left  and  right  circular  mixing  produces  linear  polarization  that  rotates  as  a  function  of 
the  transverse  focal  position.  This  rotation  is  purely  sinusoidal  and  consequently  produces  just  two  output  beams,  one  left  circularly  polarized  and  one  right 
circularly  polarized. 


polarized  HHG.  Finally,  numerical  simulations  show  that 
this  non-collinear  generation  method  extends  to  generating  an 
isolated  attosecond  burst.  To  our  knowledge,  this  is  the  first 
method  capable  of  producing  an  isolated  attosecond  burst  with 
circular  polarization. 

NCP-HHG 

To  generate  circularly  polarized  high  harmonics  in  a  non-collinear 
geometry,  two  tightly  focused,  counter-rotating  45  fs  laser  pulses 
were  overlapped  in  a  gas  jet  of  noble  gas  (xenon,  argon  or  neon), 
where  the  HHG  process  took  place  (Methods).  The  generated 
EUV  beams  propagated  in  vacuum  and  passed  through  an  alu¬ 
minium  filter  (200-nm-thick,  20-72  eV  transmission  window), 
which  blocked  any  scattered  visible  light.  The  EUV  beams  were 
then  dispersed  by  an  imaging  spectrometer  and  recorded  with  an 
X-ray  charge- coupled  device  (CCD)  camera.  Harmonics  were  gen¬ 
erated  only  when  the  beams  were  overlapped  both  temporally  and 
spatially.  Otherwise,  the  HHG  mechanism  was  suppressed  due  to 
the  circular  polarization  of  the  individual  laser  beams.  Two  EUV 
beams  were  produced  for  each  odd  harmonic  order,  with  one  circu¬ 
larly  polarized  beam  propagating  to  the  right  of  the  centreline  and 
the  other  (with  opposite  helicity)  to  the  left  (Fig.  2a). 

When  two  800  nm  driving  beams  were  used  to  generate  harmo¬ 
nics  in  argon  gas,  the  resulting  high  harmonic  photon  energies 
reached  the  29th  harmonic  order  (45  eV).  In  xenon  gas,  only 
100  |ij  in  each  driving  laser  beam  was  sufficient  to  generate  high 
harmonics,  demonstrating  that  the  NCP-HHG  method  is  easily 
scalable  to  current  multi-kHz  laser  systems.  Numerical  simulations 
(Methods)  confirmed,  for  the  case  of  perfect  alignment,  that  the  har¬ 
monics  are  fully  circularly  polarized  with  opposite  helicity  (Fig.  2b, c). 

The  nature  of  the  harmonic  emission  can  be  adjusted  by  chan¬ 
ging  the  wavelength  of  the  driving  laser.  When  two  400  nm  pulses 
were  used  to  drive  the  process  in  argon  (Supplementary  Figs  3 
and  4),  a  single  high-flux  (~2  x  108  photons  per  pulse)  harmonic 


at  22  eV  was  observed,  and  two  harmonics  (22  and  28  eV)  were 
generated  when  neon  was  used  as  the  HHG  medium.  When  267  nm 
beams  were  used  to  drive  the  NCP-HHG  process,  angularly  separated 
beams  at  14  and  23  eV  were  generated  (Fig.  la,  Supplementary  Fig.  2). 
Note  that  the  NCP-HHG  method  produces  circularly  polarized  beams 
that  are  separated  by  helicity  and  harmonic  order,  in  contrast  to 
methods  such  as  attosecond  lighthouse27,28  or  non-collinear 
optical  gating29,30,  which  produce  angularly  separated  attosecond 
pulses  with  linear  polarization  (Supplementary  Section  1). 

In  NCP-HHG,  the  emission  angle  of  the  two  EUV  HHG  beams 
is  a  small  fraction  of  the  crossing  angle  of  the  driving  laser  beams, 
resulting  in  harmonics  that  are  angularly  separated  from  the 
driving  laser.  The  divergence  angle  6q  of  the  harmonic  beams 
decreases  with  increasing  harmonic  order  q  and  closely  follows 
the  trend  sin  6q  =  ±sin  QJq,  where  Qx  is  the  crossing  half- angle  of 
the  driving  laser  beams  (Fig.  2d).  The  angular  dependence  of  the 
harmonic  emission  can  be  explained  by  a  simple  photon  model 
that  takes  into  consideration  the  conservation  of  parity,  linear 
momentum  and  angular  momentum  (Fig.  3a, c)  for  the  photons 
absorbed  and  emitted  in  the  HHG  process.  Because  of  parity 
conservation,  the  number  of  photons  absorbed  to  generate  a  given 
harmonic  n  =  n+  +  n_  must  be  odd  (n+  and  n_  refer  to  the  number  of 
photons  absorbed  from  each  of  the  two  pump  beams).  Conservation 
of  linear  momentum  dictates  that  the  output  momentum  of  the 
harmonic  photon,  k,  is  related  to  the  momentum  of  the  input 
photons  (k+  and  k_)  by  simple  vector  addition  k=  n+ k+  +  n_ k_.  As 
these  are  the  only  two  requirements  for  linearly  polarized  beams, 
the  non-collinear  mixing  of  two  linearly  polarized  beams  generates 
harmonics  at  a  variety  of  divergence  angles  (Fig.  3a),  including  both 
angles  smaller  than  the  crossing  angle  of  the  driving  beams  (sum 
frequency  generation)  and  angles  larger  than  the  crossing  angle  of 
the  driving  beams  (difference  frequency  generation)31. 

In  contrast,  in  the  case  of  counter-rotating  (left  and  right)  circu¬ 
larly  polarized  driving  beams,  conservation  of  spin  angular 

3 


NATURE  PHOTONICS  |  ADVANCE  ONLINE  PUBLICATION  |  www.nature.com/naturephotonics 


©  2015  Macmillan  Publishers  Limited.  All  rights  reserved 


ARTICLES 


NATURE  PHOTONICS  DOI:  10.1038/NPHOTON. 2015.181 


a 
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Figure  4  |  NCP-HHG  driven  by  different  frequency  driving  lasers  (400  and  800  nm).  a,  The  angle-integrated  harmonic  spectrum  reveals  the  left,  right, 
missing'  pattern  seen  in  collinear  circularly  polarized  HHG21'24.  a.u.,  arbitrary  units,  b,  In  contrast  to  collinear  generation,  non-collinear  HHG  generates  one 
helicity  in  one  direction  and  the  opposite  helicity  in  the  other  direction.  The  low  intensity  peaks  at  the  15th,  18th  and  21st  harmonic  orders  are  a  result  of  the 
slight  ellipticity  of  the  driving  laser  beams. 


momentum  adds  an  additional  constraint  on  the  allowed  combi¬ 
nations,  namely  that  nR=nL±  1.  This  selection  rule,  combined 
with  conservation  of  momentum21,32,33,  dictates  that  a  given  harmo¬ 
nic  can  only  produce  two  beams  (Fig.  3c).  The  output  beam  direc¬ 
tions  for  harmonic  order  q  follow  sin  6q  =  ±sin  OJq,  in  agreement 
with  the  experiment  and  numerical  simulations  (Fig.  2d). 
Significantly,  all  the  harmonics  are  directed  away  from  the  pump 
beams  and,  provided  the  initial  crossing  angle  is  sufficiently  large 
and  the  harmonic  divergence  sufficiently  small,  individual  harmo¬ 
nics  can  be  separated  without  using  a  diffraction  grating  (Fig.  la 
and  Supplementary  Fig.  2). 

Further  insight  into  the  non-collinear  HHG  process  can  be 
gained  by  considering  non-collinear  mixing  within  a  wave  model. 
The  crossing  of  two  linearly  polarized  beams  (Fig.  3b)  produces 
an  interference  pattern  at  the  focal  plane  and,  because  the  high  har¬ 
monic  yield  is  highly  nonlinear  with  respect  to  the  driving  field, 
there  are  many  diffracted  orders  in  the  far  field.  In  contrast,  the 
crossing  of  two  circularly  polarized  counter-rotating  beams 
(Fig.  3d)  does  not  produce  an  interference  pattern  at  the  crossing 
plane  because  the  polarizations  are  orthogonal.  However,  provided 
that  the  amplitudes  of  the  beams  are  equal,  the  combined  polariz¬ 
ation  is  locally  linear.  Because  the  relative  phase  of  the  beams 
varies  across  the  focal  spot,  the  orientation  of  the  linear  polarization 
also  varies  (Figs  lc  and  3d): 

E(x)  =  E0e~x  lw  [  cos(/c1x  sin  Q1 )  -  sin(/c1x  sin  61 )  ] 

where  x  is  the  spatial  coordinate  across  the  focal  spot  and  co  deter¬ 
mines  the  beam  waist.  The  linear  polarization  that  rotates  across  the 
focal  spot  forms  a  ‘rotating  polarization  grating  (Figs  lc  and  3d), 
which  is  the  key  to  NCP-HHG. 

In  the  interaction  region,  the  atoms  emit  high  harmonic  radi¬ 
ation  polarized  in  the  direction  of  the  local  (linearly  polarized)  fun¬ 
damental  field.  We  consider  the  propagation  of  the  horizontal  and 
vertical  components  of  the  source  polarization  separately.  The  hori¬ 
zontally  polarized  component  can  be  treated  as  a  plane  wave  that  is 
phase  modulated  with  a  spatial  cosine  pattern.  In  the  far  field,  there 
will  be  only  two  diffracted  orders  (+1,-1),  with  no  zero-order  (d.c.) 
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component  (Fig.  3d,  top).  Similarly,  the  vertically  polarized  com¬ 
ponent  is  also  phase  modulated,  but  with  a  sine  spatial  pattern. 
From  Fourier  optics,  the  90°  spatial  phase  shift  of  this  emitting 
pattern  leads  to  a  ±90°  phase  shift  of  the  ±1  diffracted  orders. 
Thus,  the  combined  fields  in  the  two  diffracted  directions  will  be  cir¬ 
cularly  polarized  with  opposite  helicity.  The  angle  of  diffraction  for 
the  qth  harmonic  is  determined  by  the  modulation  period  of  the  fun¬ 
damental  and  is  given  by  sin  6q  =  sin  OJq,  in  agreement  with  the 
photon  model. 

The  locally  linear  field  in  the  co  +  co  NCP-HHG  method  differen¬ 
tiates  this  technique  from  co  +  2co  mixing  schemes21,24,  which 
produce  a  trefoil  (three-leaf  clover)  shaped  field.  Because  of  the 
two-dimensional  nature  of  the  electric  field  in  the  co  +  2co  method, 
the  field- driven  electron  wave  packet  is  driven  back  to  the  ion 
from  different  directions  as  a  function  of  time22,23,34,  and  the 
kinetic  energy  of  the  electron  at  recollision  is  lower  than  a  linear 
field  at  co.  The  mechanism  at  play  in  the  co  +  co  non-collinear 
scheme  is  qualitatively  distinct,  as  each  atom  experiences  a  linearly 
polarized  driving  field.  This  is  advantageous,  because  the  same 
single- atom  physics  of  linear  HHG  apply,  allowing  high  photon 
energies  to  be  produced. 

An  additional  benefit  of  the  NCP-HHG  method  is  that  the  gen¬ 
erated  harmonics  are  angularly  separated  from  the  driving  laser 
(Fig.  la).  This  separation  does  not  occur  for  single-beam  or  collinear 
HHG,  where  the  harmonics  co-propagate  with  the  driving  laser(s) 
after  the  interaction  region.  In  this  case,  optically  dense  (typically 
aluminium)  filters  must  be  placed  between  the  HHG  source  and 
the  sample  to  block  the  driving  laser  and  prevent  sample  damage. 
The  filters  can  attenuate  much  of  the  EUV  (often  by  a  factor  of 
20  or  more)  and  require  that  the  sample  be  placed  far  from  the  inter¬ 
action  region  (typically  40  cm  or  more),  allowing  the  harmonic 
beam  to  diverge  significantly.  In  contrast,  using  NCP-HHG,  we 
demonstrate  that  a  fragile  sample  (in  this  case,  a  single  human 
hair)  can  be  placed  just  1  cm  from  the  interaction  region  without 
damage  (Supplementary  Fig.  4d).  In  this  geometry,  the  sample 
experiences  an  unusually  high  EUV  intensity  (on  the  order  of 
GW  cm-2)  without  the  need  for  filters,  mirrors  or  gratings. 
Considering  the  high  expense  and  low  efficiency  of  optics  in  the 
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Figure  5  |  EUV  MCD  of  an  iron  film,  a,  Spectrum  of  circularly  polarized  harmonics  generated  from  800  nm  beams  focused  in  argon  and  transmitted  through 
a  20-nm-thick  iron  sample  at  45°  incidence,  b,  MCD  asymmetry  (AMCD  =  (/+-/_)/(/+  +  /_))  between  opposing  magnetization  states  in  the  iron  sample.  As 
expected,  left  circularly  polarized  harmonics  (red)  display  opposite  dichroism  to  right  circularly  polarized  harmonics  (blue),  c,  Integrated  MCD  asymmetry 
and  95%  confidence  interval  demonstrate  that  the  dichroism  is  antisymmetric  for  the  left  and  right  polarizations  and  is  in  good  quantitative  agreement  with 
the  predicted  magneto-optical  contrast  based  on  the  available  synchrotron  data39'40  (Supplementary  Section  6). 


EUV  spectral  region,  this  simplification  of  the  apparatus  constitutes 
a  substantial  advantage. 

Non-collinear  HHG  phase  matching 

Experimentally,  the  NCP-HHG  method  generates  a  high  flux  of 
EUV  photons,  which  is  remarkable  considering  that  the  length  of 
the  interaction  region  is  limited  when  compared  to  single-beam 
or  collinear  HHG.  Fortunately,  the  decreased  length  of  the  inter¬ 
action  region  is  partially  compensated  by  a  modification  of  the 
phase-matching  conditions  that  arises  due  to  the  non-collinear  geo¬ 
metry.  Although  the  microscopic  (single-atom)  physics  of  the  circu¬ 
larly  polarized  non-collinear  HHG  method  is  identical  to  that  for 
single-beam  linearly  polarized  HHG,  macroscopically  the  crossed- 
beam  mixing  adds  a  geometric  component  to  the  phase-matching 
equation.  Specifically,  for  non-collinear  mixing  of  plane  waves,  we 
take  the  projection  of  the  k- vector  of  the  pump  beams  along  the 
direction  of  the  output  harmonic  beam: 

A k  =  nLk1cos(01  -  6q)  +  nRk1cos(01  +  6q)  -  kq 

where  each  k  includes  the  bound  and  free  electron  contributions 
to  the  refractive  index  of  the  medium  at  the  corresponding  fre¬ 
quency.  The  effective  refractive  index  of  the  harmonic  light  is 
close  to  unity,  but  the  refractive  index  at  the  fundamental  is 
greater  than  the  one  for  bound  electrons  and  decreases  in  the 


presence  of  free  electrons.  Provided  that  the  free  electron  density 
is  below  a  critical  value3,35,  the  target  density  can  be  increased  to 
achieve  phase  matching,  just  as  it  can  for  single-beam  HHG.  This 
situation  holds  for  any  crossing  angle,  because  the  bound  and  free 
electron  contributions  to  the  refractive  index  are  both  multiplied 
by  the  cosine  angular  factor. 

Thus,  as  a  result  of  phase  matching,  the  optimal  pressure  for  har¬ 
monic  generation  increases  dramatically  with  increasing  crossing 
angle.  This  increase  in  pressure  tends  to  compensate  for  the 
decreased  length  of  the  interaction  region  that  is  imposed  by  the 
non-collinear  geometry.  As  a  result,  NCP-HHG  can  reach  high 
flux  levels  as  long  as  sufficient  gas  pressures  can  be  provided.  In 
this  initial  experiment,  the  maximum  achievable  pressure  in  the 
interaction  region  was  limited  to  ~1  kPa  by  the  capacity  of 
the  vacuum  pump,  limiting  phase  matching  to  the  EUV  range.  In 
the  future,  high-pressure  gas  targets  should  allow  NCP-HHG  to 
extend  to  soft  X-ray  wavelengths. 

Non-collinear  HHG  with  different  frequency  driving  lasers 

The  circularly  polarized  non-collinear  HHG  process  is  not  limited 
to  mixing  two  beams  of  the  same  frequency.  We  demonstrated 
this  by  generating  circularly  polarized  harmonics  by  mixing  400 
and  800  nm  beams  in  argon  (Fig.  4).  The  angle-integrated  spectrum 
(Fig.  4a)  of  the  EUV  emission  exhibits  the  familiar  ‘left,  right, 
missing5  pattern  seen  in  collinear  circularly  polarized  HHG,  where 
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Figure  6  |  Numerical  simulations  reveal  additional  capabilities  of  the  NCP-HHG  method,  a,  When  driven  by  4.3-cycle,  400  nm  pulses  at  a  large  (64  mrad) 
crossing  angle,  angularly  separated  harmonics  are  produced,  obviating  the  need  for  a  spectrometer,  b,  In  the  time  domain,  each  of  these  harmonics  manifests 
as  an  angularly  separated  pulse  with  a  timescale  similar  to  the  driving  laser,  c,  When  driven  with  1.1-cycle  800  nm  pulses  (32  mrad  crossing  angle), 
supercontinuum  harmonics  are  produced,  d,  In  the  time  domain,  the  EUV  emission  manifests  as  two  isolated  attosecond  pulses,  one  left  circularly  polarized 
and  one  right  circularly  polarized.  Colour  scales  are  in  linear  arbitrary  units.  Each  panel  is  normalized  separately. 


the  left  circular  harmonics  are  separated  in  frequency  from  the  right 
circular  harmonics,  and  every  third  harmonic  is  sup¬ 
pressed21,23,24,36’37.  Furthermore,  the  angle-resolved  spectrum 
(Fig.  4b)  reveals  that  the  left  circularly  polarized  harmonics  are 
emitted  in  one  direction,  and  the  right  circular  harmonics  in  the 
opposite  direction.  Moreover,  in  contrast  to  non-collinear  mixing 
at  the  same  frequency,  where  the  harmonics  are  grouped  around  the 
centreline  of  the  pump  beams,  the  harmonics  generated  from  the 
mixing  of  400  and  800  nm  beams  are  displaced  towards  the  direction 
of  the  400  nm  output  beam  because  of  momentum  conservation 
(Supplementary  Section  4). 

EUV  MCD  measurements 

One  of  the  most  powerful  uses  for  ultrashort  pulses  of  circularly 
polarized  EUV  and  soft  X-ray  light  is  the  study  of  femtosecond 
dynamics  of  magnetic  materials  using  EUV  and  X-ray  MCD38.  In 
these  measurements,  the  magnetic  state  of  a  material  is  probed  by 
the  dichroic  absorption  of  circularly  polarized  light,  which  can 
provide  both  elemental  specificity  and  distinction  between  the 
spin  and  orbital  contributions  to  the  atomic  magnetic  moment. 
Because  the  dichroic  absorption  scales  linearly  with  the  degree 
of  circular  polarization  (Methods),  MCD  can  also  provide  a 
rigorous  measurement  of  the  ellipticity  and  helicity  of  the  EUV 
light.  The  dichroism  can  be  measured  either  by  comparing  the 
difference  in  absorption  between  the  left  and  right  circularly 
polarized  light  or  by  using  a  magnetic  field  to  flip  the  sample 
magnetization.  We  performed  an  EUV  MCD  measurement  by 
switching  the  magnetic  field  applied  to  a  20-nm-thick  iron  film 
and  comparing  the  transmitted  EUV  intensities  1+  for  opposing 
magnetization  directions,  as  expressed  by  the  MCD  asymmetry 
AMcd  =  (W-)/U+  +  U 

Experimentally,  the  MCD  asymmetry  is  readily  apparent,  even  in 
the  two-dimensional  image  (Fig.  5b),  and  clearly  demonstrates  that 
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the  two  beams  exhibit  opposite  helicity.  Angular  integration 
(Fig.  5c)  provides  a  quantitative  estimate  of  the  MCD  asymmetry 
of  ~1.5%  in  the  range  of  29-39  eV.  To  our  knowledge,  there  is  no 
experimental  measurement  of  the  MCD  contrast  (or  magneto¬ 
optical  constants)  of  iron  in  this  energy  range,  but  synchrotron 
data  at  slightly  higher  photon  energies  are  available39,40.  The  trend 
of  the  data  indicates  that  the  MCD  asymmetry  in  the  range  of  our 
measurement  is  less  than  2%  (Supplementary  Fig.  6).  Thus,  our 
measurement  of  ~  1.5%  asymmetry  confirms  that  the  EUV  light  gen¬ 
erated  by  the  NCP-HHG  process  has  a  high  degree  of  circular  polar¬ 
ization.  Furthermore,  this  EUV  MCD  measurement  demonstrates 
that  the  NCP-HHG  source  is  a  practical  tool  for  measuring  ultrafast 
magnetic  dynamics.  The  data  presented  in  Fig.  5  represent  30  min  of 
averaging,  but  an  adequate  measure  of  MCD  asymmetry  can  be 
acquired  in  2  min  (Supplementary  Fig.  7),  opening  the  door  to 
time-resolved  pump-probe  studies  of  femtosecond- resolved 
magnetic  dynamics. 

Future  capabilities  investigated  with  numerical  simulations 

Numerical  simulations  of  HHG  including  propagation  (Methods) 
provide  insight  into  further  capabilities  of  the  NCP-HHG 
method.  First,  experimentally,  the  harmonics  could  only  be  angu¬ 
larly  separated  when  using  267  nm  driving  lasers  due  to  the  low 
pressure  in  the  interaction  region  (Fig.  1).  However,  simulations 
predict  that  similar  separation  can  be  achieved  at  longer  driving 
laser  wavelengths  and  higher  harmonic  photon  energies.  For 
example,  simulations  indicate  that  400  +  400  nm- driven  NCP- 
HHG  at  a  64  mrad  crossing  angle  will  produce  harmonics  that  are 
completely  separated  angularly  (Fig.  6a).  The  situation  is  similar 
for  longer  wavelengths,  but  the  required  crossing  angle  becomes 
larger.  The  separation  of  the  harmonics  allows  for  spectroscopic 
experiments  such  as  X-ray  transient  absorption  spectroscopy  and 
MCD  to  be  completed  without  a  spectrometer  to  disperse  the 
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harmonics,  a  distinct  advantage  considering  the  high  expense,  low 
efficiency  and  temporal  dispersion  of  gratings  in  the  EUV  and 
soft  X-ray  regions.  In  this  case,  the  HHG  process  itself  serves  as 
the  dispersive  grating  (Fig.  3d).  Of  course,  where  separation  of  the 
harmonic  beams  is  not  desirable,  a  small  crossing  angle  can  be 
used  to  achieve  good  overlap  of  the  harmonic  orders,  and  the 
NCP-HHG  method  can  still  serve  as  a  practical  source  of  circularly 
polarized  high  harmonic  light. 

Simulations  also  provide  clear  insight  into  the  temporal  structure 
of  the  EUV  emission.  When  all  the  harmonics  are  angularly  separ¬ 
ated  (Fig.  6a),  the  EUV  beams  emerge  as  single  long  pulses  (Fig.  6b), 
with  a  temporal  duration  somewhat  shorter  than  that  of  the  driving 
laser  (that  is,  a  several-femtosecond  pulse).  When  the  harmonics  are 
angularly  overlapped,  the  standard  ‘attosecond  pulse  train’  is  gener¬ 
ated,  which  consists  of  numerous  attosecond  bursts  (Fig.  2e). 
However,  with  few- cycle  driving  pulses,  angularly  dispersed  super¬ 
continuum  harmonics  are  generated  (Fig.  6c).  In  the  time  domain, 
this  manifests  as  an  isolated  attosecond  pulse  of  circularly  polarized 
EUV  light  (Fig.  6d).  Though  previous  theoretical  works  have  pro¬ 
posed  the  generation  of  isolated  circularly  polarized  attosecond 
pulses  using  exotic  conditions,  such  as  very  intense  terahertz 
fields41,  to  our  knowledge,  the  NCP-HHG  method  is  the  first 
experimentally  realized  method  that  is  capable  of  producing  isolated 
attosecond  pulses  with  circular  polarization. 

Discussion 

Looking  forward,  the  NCP-HHG  method  is  poised  to  enable  break¬ 
through  probes  of  ultrafast  dynamics.  When  driven  with  few-cycle 
pulses,  simulations  indicate  that  the  NCP-HHG  method  will 
produce  two  circularly  polarized  isolated  attosecond  pulses,  one 
propagating  to  the  left  and  one  to  the  right.  This  geometry  naturally 
lends  itself  to  an  attosecond  pump-attosecond  probe  experiment, 
which  could  be  implemented  using  a  single  ‘split-mirror’  focusing 
optic  to  adjust  the  temporal  delay.  Because  the  NCP-HHG 
method  can  eliminate  the  need  for  all  other  optics,  gratings  and 
filters,  both  the  EUV  intensity  at  the  sample  and  the  EUV  flux  at 
the  detector  could  reach  levels  unprecedented  in  attosecond 
pump-attosecond  probe  experiments,  potentially  transforming 
these  formidable  experiments  into  commonplace  techniques  for 
probing  the  fastest  processes  in  materials,  molecules  and  atoms. 
Similarly,  for  imaging  studies,  refocusing  the  two  beams  of  opposite 
helicity  with  a  single  optic  will  allow  for  new  control  over  illumina¬ 
tion,  proving  either  separated  right  and  left  circular  polarization,  or 
a  ‘textured’  linear  polarization  reproducing  the  rotating  polarization 
pattern  of  the  source.  Furthermore,  recent  work  has  shown  that 
high-intensity  circularly  polarized  EUV  pulses  can  be  used  to  generate 
attosecond  magnetic  field  pulses42. 

Here,  we  have  demonstrated  NCP-HHG  for  the  first  time, 
and  have  shown  that  this  method  generates  bright  circularly 
polarized  EUV  beams  of  both  left  and  right  helicity  simultaneously. 
Owing  to  selection  rules,  the  two  helicities  are  angularly  isolated, 
and  each  harmonic  order  is  emitted  at  a  different  angle  (Fig.  1), 
allowing  the  HHG  process  itself  to  serve  as  a  spectrometer. 
Conveniently,  the  EUV  light  is  well  separated  from  the 
intense  driving-laser  beams.  Using  this  new  light  source,  we 
performed  an  EUV  MCD  measurement  that  confirmed  the  circular 
polarization  of  the  EUV  light  and  demonstrated  that  the  flux  and 
stability  of  the  source  are  sufficient  for  real-world  studies  of 
femtosecond  magnetic  dynamics.  Importantly,  we  demonstrated 
that  the  single-atom  physics  of  the  NCP-HHG  process  is  identical 
to  that  of  single-beam  linearly  polarized  HHG,  indicating  that  the 
NCP-HHG  method  can  be  extended  to  produce  circularly  polarized 
harmonics  in  the  soft  X-ray  region  by  using  mid-infrared  driving 
lasers43,44,  produce  high-flux  EUV  sources  using  ultraviolet 
driving  lasers45,  and  generate  isolated  attosecond  pulses  using 
few-cycle  driving  lasers46. 


Methods 

Methods  and  any  associated  references  are  available  in  the  online 
version  of  the  paper. 
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Methods 

Non-collinear  HHG  experiment.  The  apparatus  for  the  generation  of 
circularly  polarized  HHG  light  consisted  of  a  Mach-Zehnder  interferometer,  which 
focused  two  laser  pulses  into  a  gas  jet  in  a  non-collinear  geometry.  The  laser  was  a 
Ti:sapphire  regenerative  amplifier  (KMLabs  Wyvern-HE)  operating  at  1  kHz  with  a 
45  fs  pulse  duration  and  a  central  wavelength  of  795  nm  (referred  to  as  ‘800  nm’). 
A  pulse  energy  of  200-400  pj  in  each  arm  was  typically  used  for  HHG  in  argon  gas. 
The  laser  pulses  were  focused  with  the  same  lens  (/=  250  mm)  into  the  gas  jet. 

Circularly  polarized  HHG  was  demonstrated  using  four  wavelength 
combinations:  800  +  800,  800  +  400,  400  +  400  and  267  +  267  nm.  The  second 
harmonic  of  the  fundamental  (400  nm)  was  generated  using  a  200-pm-thick  beta 
barium  borate  (BBO)  crystal.  The  third  harmonic  of  the  fundamental  (267  nm)  was 
generated  by  four-wave  mixing  of  the  400  and  800  nm  beams  in  a  second 
BBO  crystal  (Eksma  Optics).  Half  and  quarter  achromatic  waveplates  placed  in  each 
arm  of  the  interferometer  were  used  to  manipulate  the  polarization  and  produce  left 
and  right  circularly  polarized  laser  pulses.  The  peak  intensity  in  the  overlapping 
beams  was  similar  to  that  required  in  single  beam  HHG,  ~1  x  1014  W  cm-2  for  argon. 
The  two  non-collinear  beams  have  a  half-angle  separation  of  1.25-2.5°  (25-50  mrad) 
and  were  overlapped  over  a  150-pm-diameter  gas  jet  nozzle.  The  HHG  media 
were  argon  and  neon  gas,  with  backing  pressures  between  100  and  300  kPa.  The 
generated  EUV  light  propagated  in  vacuum  (<0.1  Pa)  and  passed  through  a  200  nm 
aluminium  filter  (Luxel)  before  being  dispersed  using  a  grating-based  spectrograph 
(Hettrick  Scientific)  and  recorded  using  a  back-illuminated  X-ray  CCD  camera 
(Andor  Newton  DO920N-BN-9HC).  The  EUV  emission  was  easily  observable 
even  with  short  exposure  times  (the  longest  exposure  time  used  in  this  study  was 
only  10  s). 

EUV  MCD  measurement.  The  sample  for  the  EUV  MCD  experiment  was 
prepared  by  first  sputtering  5  nm  germanium  onto  a  500-nm-thick 
unsupported  aluminium  filter.  On  top  of  the  germanium,  20  nm  untextured  iron 
was  deposited,  which  formed  the  actual  magnetic  sample  to  be  measured.  The 
iron  layer  was  capped  with  another  5  nm  germanium  to  prevent  oxidation. 

For  the  MCD  measurement,  an  electromagnet  was  used  to  apply  a  15  mT  magnetic 
field  for  saturating  and  switching  the  magnetization  of  the  sample 
(Supplementary  Fig.  5).  After  collecting  a  10  s  exposure,  the  sample  magnetization 
was  switched  and  the  data  collection  was  repeated  for  a  total  of  90  exposures  in  each 
field  direction. 

The  resonant  magneto-optical  effect  can  be  described  as  helicity-dependent 
absorption  in  a  magnetic  material,  where  the  imaginary  (absorptive)  part  of  the 
refractive  index  n+  of  the  magnetically  saturated  material  is  given  by  Im(n±)  =  />  ±  A/1, 
with  />  denoting  the  average  absorption  coefficient  and  A/1  ~  Pcircmcos  a  the  dichroic 
contribution  for  opposing  photon  helicities  (±),  where  Pcirc  is  the  degree  of  circular 
polarization  of  the  light  and  me  os  a  is  the  projection  of  the  magnetic  moment  in  the 
direction  of  the  incident  light47,48.  Thus,  the  MCD  asymmetry  can  be  measured  by 
switching  the  magnetic  field  with  an  external  electromagnet.  To  this  end,  a  20-nm- 
thick  iron  sample  was  inserted  at  an  angle  of  incidence  a  of  45°  into  the  HHG  beams 
(Supplementary  Section  6).  The  in-plane  magnetization  of  the  sample  was  switched 
with  an  electromagnet  and  the  integrated  signal  from  opposing  sample 
magnetizations  was  recorded  with  an  X-ray  CCD  camera. 

Theoretical  HHG  simulations  including  propagation.  The  harmonic  propagation 
was  computed  using  a  method  based  on  the  electromagnetic  field  propagator49.  The 
target  (gas  jet)  was  discretized  into  elementary  radiators  and  propagated  the  emitted 


field  E;(rd,f)  to  the  detector, 


where  qj  is  the  charge  of  an  electron,  sd  is  the  unitary  vector  pointing  to  the  detector, 
and  rd  and  r;  are  the  position  vectors  of  the  detector  and  the  elementary  radiator  j, 
respectively.  The  dipole  acceleration  a;  of  each  elementary  source  was  computed 
using  the  SFA+  method50,  an  extension  of  the  strong  field  approximation51-53.  The 
signal  at  the  detector  was  computed  as  the  coherent  addition  of  the  HHG 
contributions  of  all  elementary  sources,  with  the  HHG  light  assumed  to  propagate  to 
the  detector  with  phase  velocity  c.  The  propagation  effects  of  the  fundamental  field, 
including  plasma  and  neutral  dispersion  as  well  as  time-dependent  group  velocity 
walk-off,  were  all  taken  into  account.  The  absorption  of  the  harmonics  in  the  gas  was 
modelled  using  Beer’s  law.  The  attosecond  pulses  were  computed  by  Fourier 
transformation  of  the  harmonic  spectrum  detected  at  a  particular  angle  after  filtering 
out  the  low-energy  portion  (lower  than  the  ninth  harmonic  order  of  the  fundamental). 

Simulations  were  performed  for  800  +  800  and  400  +  400  nm  counter-rotating 
non-collinear  mixing,  with  32  and  64  mrad  crossing  angles,  respectively.  The  laser 
pulse  envelopes  were  modelled  by  a  sin2  function  truncated  at  the  first  zeros,  with  a 
peak  intensity  of  8.75  x  1013  W  cm-2  at  the  focus  for  each  beam.  Various  pulse 
durations  were  simulated  to  investigate  the  role  of  pulse  duration  in  the  HHG 
process.  Each  non-collinear  driving  beam  was  assumed  to  have  a  Gaussian  spatial 
profile  (with  a  beam  waist  of  40  pm)  and  taken  directly  as  an  input  to  the 
calculations.  The  target  was  modelled  as  an  argon  gas  jet  flowing  perpendicular  to 
the  longitudinal  direction  z.  The  gas  jet  profile  was  constant  in  the  x  direction,  but 
Gaussian  in  the  y  and  z  directions,  with  full-width  at  half-maximum  (FWHM) 
values  of  40  and  20  pm,  respectively,  and  a  peak  pressure  of  667  Pa.  Note  that  the 
width  of  the  gas  jet  along  the  z  direction  was  sufficiently  small  to  assume  linear 
polarization  of  the  total  electric  field  along  the  target. 

Note.  The  suppliers  of  commercial  equipment,  instruments  or  materials  are 
identified  in  this  Article  to  foster  understanding.  Such  identification  does  not  imply 
recommendation  or  endorsement  by  the  National  Institute  of  Standards  and 
Technology,  nor  does  it  imply  that  the  materials  or  equipment  identified  are 
necessarily  the  best  available  for  the  purpose. 
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